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We designed and built a new Four-Stokes-Parameter spectral line Polarimeter (FSP- 
Pol) for the Caltech Submillimeter Observatory (CSO). The simple design of FSPPol 
does not include any mirrors or optical components to redirect or re-image the radiation 
beam and simply transmits the beam to the receiver through its retarder plates. FSPPol 
is currently optimized for observation in the 200 — 260 GHz range and measures all four 
Stokes parameters, /, Q, U, and V. The very low level of instrument polarization makes 
it possible to obtain reliable measurements of the Goldreich-Kylafis effect in molecular 
spectral lines. Accordingly, we measured a polarization fraction of a few percent in 
the spectral line wings of 12 CO (J = 2 — > 1) in Orion KL/IRc2, which is consistent 
with previous observations. We also used FSPPol to study the Zeeman effect in the 
N = 2 — > 1 transition of CN in DR21(OH) for the first time. At this point we cannot 
report a Zeeman detection, but more observations are ongoing. 

Subject headings: Instrumentation: Polarimeters — ISM: Clouds — ISM: Magnetic Fields — 
Submillimeter 
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1. Introduction 

Astronomical polarimetry is a powerful tool for studying the characteristics of the interstellar 
medium, from the large scale galactic magnetic fields to the gravitational collapse of molecular cloud 
cores. The understanding of the physical phenomena responsible for polarized interstellar emission 
provides valuable information about the underlying astrophysical processes, an important example 
being the formation and evolution of stars. Since the 1990s, polarimetry tech niques and instruments 



have been developed for a wid e range of wavelengths, from optical (e.g. 
Wiktoro wicz fc Matthews! 120081) to the submillimeter and radio regimes fe.g 
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One of the important applications of astronomical polarimetry lies in the study of the effect 
of magnetic fields on the early stages of star formation. Hence, it is important to characterize the 
interstellar magnetic fields around star-forming regions. Indeed, it is possible to obtain the orienta- 
tion and strength of the magnetic fields observationally, as they leave their signature on the emission 
of the interstellar dust grains and mol ecules. For example, th e dust continuum radiation becomes 



linearly polarized in a magnetic field (jHildebrand et al. 



1999). At submillimeter wavelengths this 



polarization is perpendicular to the field lines and a polarimetry map of dust continuum emission 
will therefore reveal the plane-of-the-sky orientation of the magnetic field. Furthermore, in the 
presence of anisotropic radiation, ma gnetic fields cause the emi ssion from gas molecules to also be 
linearly polarized by a few percent (JGoldreich fe Kvlafislll98ll ). It is also possible to determine 



whether this polarization is aligned parallel or perpen dicular to the plane-of-the-sky component of 
the ambient magnetic field (jDeguchi &; Watsonlll984l ). Finally, the only way to directly measure 
the strength of interstellar magnetic fields is through Zeeman line-broadening measurements. More 
precisely, the line-of-sight component of the field is obtained by measur ing the circular polarizatio n 



in the emission of the Zeeman components of molecular spectral lines (jCrutcher et al. 



1993 



1999). 



We describe the Four-Stokes-Parameter spectral line Polarimeter (FSPPol) we recently designed 
and successfully commissioned at the Caltech Submillimeter Observatory (CSO) in November 2008 
in §2. We present linear polarization measurements of 12 CO (J = 2 — > 1) in Orion KL/IRc2 in §3 
and preliminary Zeeman observations of the N = 2 — > 1 transition of CN in DR21(OH) in §4. We 
end with a summary in §5. 



2. Instrument Description 

The design of FSPPol is based on single beam polarimetry, as the current heterodyne receivers 
at the CSO can measure only one linear polarization state at one position at a time. The simple 
construction of the polarimeter enabled us to mount the instrument inside the elevation tube that 
traces much of the optical path from the tertiary mirror behind the telescope dish to the Nasmyth 
focus where the 200 — 300 GHz receiver is located. A schematic diagram of FSPPol installed in 
the elevation tube is shown in Figure [TJ The polarimeter is mounted on a bracket attached to the 
walls of the tube and its location with respect to the tertiary mirror and the receiver is depicted 
with the schematic diagram of the CSO telescope in Figure [2] The third focus of the 230 GHz 
telescope beam is virtual and located 1.07 m behind the tertiary mirror and has a waist of 35.9 
mm. We calculated the beam waist upon incidence on FSPPol to be 39.5 mm. As illustrated 
in Figure [H FSPPol is comprised of a half- wave plate (HWP) and a quarter- wave plate (QWP), 
each being 100 mm in diameter and optimized for observations at 226 GHz. The wave plates were 
manufactured by Meller Optics, Inc^J and subsequently anti-reflection coated by QMC Instruments 
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Ltdo More detailed specifications of the wave plates are listed in Table [TJ These plates are installed 
in rotating rings mounted side by side on an aluminum translational stage that moves across the 
elevation tube, placing either wave plate in the path of the signal reflected from the tertiary mirror. 
The anti-reflection coating on the plates is efficient in eliminating standing waves along the optical 
path. Nevertheless, the mounting of the translational stage was adjusted in such a manner that 
the telescope beam impinged on the HWP or QWP at an incidence that is slightly off from normal 
(i.e., at most a few degrees) to avoid potential standing waves. The rotational and translational 
motions of these wave plates are precisely monitored and controlled by the instrument controlling 
software from the observatory's control room. The system temperature increases by about 8% when 
the HWP is placed in the path of the beam, and by about 6% with the QWP in the beam. 

We aligned the optical axes of the wave plates in relation with the polarization axis of the 
receiver by measuring the intensity from a cold load, located at the entrance of FSPPol, while 
rotating the wave plates over a wide range of angles. We previously placed a wire grid with its 
polarizing state parallel to the receiver's axis between the wave plates and the cold load to maintain 
a specific incoming polarization state for testing purposes; this grid was removed for astronomical 
observations. It is straightforward to show that the location of a minimum in measured intensity 
with this set-up corresponds to an angle where either the fast or slow axis of a plate (HWP or 
QWP) is aligned with the receiver axis. This allows for a precise orientation of the HWP and QWP 
with the receiver axis (also see §3). It should also be noted that the slow axis is marked by the 
manufacturer on each plate, which we used as a further consistency check. A small circular cold 
load was also used to position the polarimeter within the elevation tube such that the telescope 
beam was centered on the surface of the HWP and QWP. 

To test the instrument's accuracy, we observed the linear polarization in the J = 2 — > 1 
transition of 12 CO towards the Orion KL/IRc2 high mass star forming region. After successful linear 
polarization tests, we proceeded to perform Zeeman observations on the JV = 24l transition of 
CN in DR21(OH) for the first time. The details of these observations are explained in the following 
sections. 



3. Linear Polarization Measurements 

We observed 12 CO (J = 2 ->• 1) at 230.54 GHz in three differe n t loca tions in Orion KL/IRc2 



coincident with previous polarization observations by iGirart et al,l (120041 ) , in order to assess the 
performance of FSPPol. These locations, defined in the equatorial coordinate system, were offset 
(20", 20"), (20", -20") and (-20", -20") from the center of IRc2 (a = 05 h 35 m 14.50 s and <5 = 
— 05°22 30.4 , J2000.0). The observations were performed between September 27 th and October 6 th 
2009 using the 200 — 300 GHz receiver at the CSO, and the FFTS spectrometer with a bandwidth 
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of 500 MHz and a channel width of 61 kHz that corresponds to a velocity resolution of about 0.08 
km s _1 . The telescope efficiency was determined with scans on Jupiter and calculated to be ~ 65 
% for a beam width of ~ 32"(FWHM). During these observations the typical optical depth at 225 
GHz, as obtained with the CSO radiometer, was T225 ~ 0.06 and the typical system temperature 
with FSPPol in use was T sys ~ 370 K. 

Although the refractive indices for the ordinary and extraordinary rays vary significantly with 
frequency, th e birefringence of crystal Quartz appears to be relatively unchanged over a wide band 



(|Afsarlll987l ) even though there is some uncertainty in th e literature over its value (0.048 for the 



frequencies we are concerned with; iMarrone and Raoll2008l ). However, since the HWP was operated 



at approximately 230.5 GHz, i.e., more than 4 GHz away from its design frequency, we characterized 
its performance at the aforementioned frequency by measuring the power from a cold load at different 
orientations of the HWP slow axis with respect to the polarization axis of the facility receiver at 
the CSO. Figure [3] shows the plot of the power at the receiver against different angles between the 
HWP and receiver axis. The power is measured in arbitrary units and the angles are displayed 
in degrees. A polarizing grid, with its polarization state parallel to the receiver axis, was placed 
between the cold load and the HWP to maintain a specific incoming polarization state. As seen 
in the figure the HWP response is yery good at that frequency, although the power at successive 



minima is observed to slightly vary (jSavini et al.ll2009l ). Our model for the fit to the data in Figure 
[3] is a cosine function with the form 

/( ^(C^ + <CL-HL) cos(4(s _ 5)) _ (1) 

where 1(9) is the power measured by the receiver at angle 9 and 5 determines the offset of the HWP 
slow axis from the orientation defining the polarization state of the receiver. CL is the power from 
the cold load and HL is the power due to emission from the receiver that reflects back on the wires 
of the polarizing grid; The HL signal incident on the HWP is therefore polarized perpendicular 
to the CL signal. It should also be noted that there is an equal contribution to CL and HL from 
receiver noise, although we cannot quantify it with these measurements alone. The fit values for 
CL, HL and 5 are 722.2 ± 2.1, 1373.0 ± 2.9 and 1.45° ± 0.11°, respectively. 

The polarization efficiency of FSPPol for linear polarization measurements was previously de- 
termined in November 2008 by taking scans on Saturn. We placed a wire grid at the entrance of 
FSPPol with its polarization state parallel to the receiver's axis to force a precise incoming polariza- 
tion state, and observed Saturn for HWP angles 9 = 0°, 90°, 45° and 135°, where 9 is defined as the 
angle between the HWP slow axis and the receiver polarization axis. We performed a one-minute 
(on-source) integration for each measurement, with the system temperature being calibrated before 
every scan and the signal from Saturn being integrated over the whole available spectral bandwidth. 
Assuming that the wire grid is perfectly polarizing and that Saturn is unpolarized, we obtained a 
polarization efficiency of ~ 99%. These measurements also allowed us to estimate the instrument 
polarization to be on the order of ~ 0.3%. In view of this high polarization efficiency, we did not 



correct our Orion KL/IRc2 data for the efficiency and instrumental effects. 

The polarimeter and the receiver are mounted in such a way that they rotate in elevation with 
the telescope and the polarization axis of the receiver on the sky, which is precisely oriented east- 
west (parallel to the horizon) when FSPPol is not used, is preserved regardless of the pointing of 
the telescope since the CSO telescope has an alt-azimuth mount. Our aim for linear polarization 
measurements was to obtain the Stokes Q and U parameters in the reference frame of the sky, i.e., 
in the equatorial coordinate system, as well as the Stokes /. Therefore, the rotation of the object's 
orientation on the sky, or in other words the changes in its parallactic angle with time had to be 
considered. For this purpose, we obtained one-minute (on-source) long intensity measurements for 
different orientations of the HWP slow axis relative to the receiver polarization axis 8 = (7 + 8 + 
90°)/2, where 7 is the parallactic angle defining the orientation of the object on the sky and 8 is the 
angle at which we seek to measure the linear polarization state in the reference frame of the source 
(i.e., the frame that is rotated by 7 in the equatorial system; see Figure U]). These measurements are 
denoted by I g i . The instrument controlling software continuously obtained the updated value of the 
parallactic angle from the observatory's antenna computer during the on-source integration, and 
compared it to the initial value at the beginning of the integration. Once the change in 7 exceeded 
a predetermined threshold (i.e., 1°), the software commanded the antenna computer to stop the 
integration, rotated the HWP by the updated angle using the new value for 7, and the integration 
resumed. Our observations were performed in cycles of four measurements for 8' = 0°, 90°, 45° and 
135°. Similar to the observations on Saturn, we took one minute long scans for each measurement, 
and calibrated the system temperature before every scan. We adopted this conservative observing 
plan to minimize the effect of pointing or calibration errors in the polarization data. The telescope 
pointing was verified hourly with typical "five points" integrations on CRL865 (a = 06 03 m 59.8 s 
and 5 = 07°25 51.4 , J2000.0) as our reference star. 

With the aforementioned definition for I a i we have 
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(Io° + -^90° + ^45° + ^135° ) 



Q — Io° — ho° (2) 

U = ^135° — ^45° > 



and 



P = j (3) 

1 (u 

PA = - arctan I — 
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for the polarization fraction and angle (measured from north increasing eastwards), respectively. 



The results from our measurements are shown in Figure [5l where the top graphs show the Stokes 
I spectra (corrected for the beam efficiency), while the middle and bottom graphs are for the 
corresponding polarization fractions and angles, respectively. For p and PA the data were binned 
using six adjacent velocity channels and only values for which p > 3a p are plotted, where a p is the 
uncertainty in the polarization fraction. 

Table [2] displays the average values for p and PA in the blue and red wings of the CO spectra, as 
well as the small polarization fraction detected at the center of the lines. These values were obtained 
after calculating the average of the Stokes /, Q and U across the stated velocity ranges. The small 
polarization level in the center of the lines is probably the contribution from the instrumental 
polarization, as it is expected that the linear polari zation due to the Gqldreic h-Kylafis effect will 

The values of PA 



be greatly reduced where the optical depth is high (IGoldreich fc Ky lafis 



1981 



are generally uniform across the spectral lines except for the center of the lines, where the relative 
contribution from the instrumental polarization is significant. 

One potential problem in the linear polarization measur ements is the presen ce of polarized 
sidelobes that introduce false polarization signals in the data (JForbrich et al.ll2008l ). This is more 
likely to happen when the telescope is pointed towards the edge of an extended source for the 
purpose of polarimetry measurements in low intensity regions. This way, it is possible that strong 
emission from the core of the source falls on these polarized sidelobes and contaminates the data. 
We have yet to determine the potential contribution of polarized sidelobes to data obtained with 
FSPPol at the CSO, we cannot therefore comment on their significance for our results on Orion KL. 
We hope to do so during a future observing run. 



Our polarization results are in general consistent with the findings of lGirart et al.l (|2004l ) (see 
their Figure 1). The polarization fraction that we calculate in the spectral line wings of 12 CO is in 
good agreement with their results, and the difference between their values for polarization angles 
and ours (i.e., 10° — 20°) could solely be the result of different telescope beam sizes. This could also 
explain the slight differences in line profiles and intensities between their Stokes / spec tra and ours. 
For ex ample the secondary peak at ~ 12 km s _1 in the offset positions obtained by iGirart et al. 
(|2004l ) is not resolved in our spectral data. 



4. Zeeman Measurements of CN (N = 2 — >■ 1) 

As mentioned earlier, measurement of the Zeeman effect in interstellar molecular spectral lines 
is the only direct way to obtain the strength of magnetic fields in molecular clouds. The Zeeman 
line-broadening in an interstellar spectral line profile is directly proportional to the strength of 
the ambient magnetic field, and can be studied by observing the signature of circular polarization 
in the line profile. To this date, Zeeman detections have been reported in the spectral lines of a 



few species, namely H I, OH, H?0, CH3OH and CN dSarma &: Momiian 



2009; 



Falgarone et al.ll2008l : ISarma et al.ll200ll : ICrutcher et al.lll999J : IPlante et al.lll995l ) . The high critical 



Troland et al. 2008 



density of CN makes this molecule suitable for studying magnetic fields in dense regions. Addition- 
ally, the rotational transition lines of CN contain several hyperfine components that have different 
Zeeman splitting coefficients. This makes it possible to distinguish between the tru e Zeeman effect 



and instrumental effects that produce artificial circular polarization in the data (jCrutcher et al. 
19961 ) . A lthough the Zeeman effect in the N = 1 — > transition of CN has previously been 



detected (IFalgarone et al.l 120081 : ICrutcher et al.l Il999l ) , to the best of our knowledge it has never 
been attempted at N = 2 —?■ 1. The higher critical density of the latter will make it possible to 
probe denser regions in molecular clouds and better establish observationally how the magnetic field 
strength scales with density. 

There are nine strong hyperfine components in the ^ = 2^1 transition of CN. The frequen- 
cies, Zeeman coefficients and relative intensities of these lines are displayed in Table [3j Similar to 
the N = 1 — > transition previously observed in DR21(OH), the CN (N = 2 — > 1) lines are double 
peaked, suggesting the existence of two velocity components (JCrutcher et al.lll999i ). These two ve- 
locity components may be arising from different regions, and may exhibit different field strengths. 
It is possible to observe the nine hyperfine components simultaneously and f it the observed circul ar 
polarization data to the following expression for both velocity components (JCrutcher et al.lll999i ) 



Vi = a(I tl + I l2 ) + biA + b 2 A + Cl (zA + c 2 (zA, (4) 

dv dv dv dv 

where Vi is the total circular polarization for hyperfine line i. In the above equation, the different 

sources of polarization are expressed in separate terms. The first term represents the polarization 

contribution from the error in intensity calibrations in the two polarization modes, a, that adds 

a small image of the intensity spectrum, In + Ii 2 , to the Vi spectra. The second and third terms 

are polarization contributions from the beam squint effect represented by b\ and b 2 , that is caused 

when the right and left handed polarization measureme nts do not probe exac tly the same region, for 

which a non-zero velocity gradient exists in its vicinity (JCrutcher et al.lll993l ) . This effect is brought 

about by pointing errors of the telescope, and also in some cases by mechanical deformations in 

the antenna. The last two terms are the true Zeeman signals for each velocity component of the 

spectrum with c\ and c 2 representing the strength of the line-of-sight component of the magnetic 

fields B\ os \/2 and B\ os2 /2, respectively, for each velocity component, and Zi being the Zeeman 

coefficient for hyperfine line i. This way, the real Zeeman effect is separated from instrumental 

effects by determining the above fitting parameters. 



4.1. Observations 

We observed CN (N = 2 ->■ 1) at 226.8 GHz in DR21(OH) (a = 20 h 39 m 01 s and 5 = 
42°22'37.7", J2000.0) during July 8 th to 13 th and September 27 th to October 6 th 2009, using the 
FFTS spectrometer with bandwidth of 500 MHz and channel resolution of 61 kHz. The telescope 
efficiency was determined with scans on Saturn and Jupiter and calculated to be ~ 65 % for a beam 



width of ~ 32"(FWHM). We verified the telescope pointing on an approximately hourly rate by 
performing 12 CO (J = 2 -> 1) scans on X Cygni (a = 19 h 50 m 33.8 s and 5 = 32°54'53.2", J2000.0) 
as our reference star. Our first observing session suffered from issues with an unstable receiver and 
mediocre skies. The average system temperature at that time was about 500 K and the typical 
r 225 was ~ 0.12. The circular polarization observations were performed by taking scans of CN 
(AT = 2 -)• 1) in DR21(OH) with the QWP slow axis rotated by +45° and -45° with respect to the 
receiver's polarization axis. With our instrument set-up, our one minute on-source integrations at 
+45° project the right-handed circular polarization (Ircp) emission (on the sky) on the receiver's 
axis, while similar integrations at —45° probe the corresponding left-handed circular polarization 
(Ilcp) emission. The Stokes / and V are obtained from I = Ilcp+Ircp an d V = Ilcp — Ircp- As 
before, this observing strategy was chosen in order to mitigate any potential calibration or pointing 
errors. The system temperature was calibrated before every scan. We have so far obtained a total 
on-source integration time of 553 minutes, with an average system temperature of 450 K. 

Figure [6] shows the Stokes I spectrum for the CN [N = 2 — > 1) hyperfine lines that are labeled 
according to the order given in Table [3] with the line temperatures corrected for the telescope beam 
efficiency. There are two further hyperfine lines at ~ —45 km s _1 and ~ 200 km s _1 that are weak 
relative to the labeled lines and therefore not included in the study. Observing the N = 2 —?■ 1 
transition of CN brings more complications compared to the N = 1 — > transition. For example, 
there are three lines in Table [3] that are labeled 5 because they are heavily blended and appear as 
a single broad line in the observed spectrum, as seen in Figure [6J Furthermore, due to the double 
peaked feature of the CN lines, the other hyperfine lines with small frequency separations (i.e., lines 
2 and 3 and lines 6 and 7) are blended together as well, although to a lesser level. 

An inspection of the ratio of line temperatures reveals that our CN lines are not consistent 
with the LTE assumption. The spacing between the two velocity components of every hyperfine 
line varies with the line strength, and the rotational diagram for the N = 2 —?■ 1 lines produced a 
negative excitation temperature, suggesting that the lines are affected by self absorption. A more 
complete spectral line analysis will be presented in a forth-coming paper. 



4.2. Preliminary Analysis & Discussion 

The following is only a preliminary analysis on the data we have gathered so far, as more data 
will be obtained in upcoming observing runs, and a more thorough discussion will be presented 
then. 

Since the CN hyperfine lines in DR21(OH) are double peaked, we fitted two to three Gaussian 
profiles to each line depending on its shape. In the fitting procedure, the line widths of the Gaussian 
profiles within each velocity component of a hyperfine line profile were kept the same. For the 
three blended hyperfine components of line 5 (in Figure [6]), we fitted one Gaussian profile to each 
component and fixed their relative frequencies according to the information given in Table [3l Each 



hyperfine component of the other pairs of blended lines (lines 2 and 3 and lines 6 and 7) were treated 
as single lines, such that Gaussian profiles incorporating both lines were fitted to the blended line 
profiles. Once the fitting parameters for all the Gaussian profiles were obtained, the blended lines 
could thus be separated by subtracting the Gaussian fit of one line from the spectrum of the pair. 
The Gaussian fits for the Stokes Jj = In + I^ 2 f° r both velocity components of every hyperfine line 
i, together with their derivatives with respect to frequency and corresponding Zeeman coefficients 
were simultaneously fitted to the Stokes V spectrum using Equation (j4]). The resulting fitting 
parameters were a = -0.0015 ± 0.0007, b\ = 1.24 ± 1.13 kHz, b 2 = 1.28 ± 1.47 kHz, c\ = -0.5 ± 1.3 
mG and c 2 = —2.5 ±2.0 mG. Although we have not obtained a Zeeman detection, the instrumental 
polarization contribution is seen to be very small for a and kept to reasonable levels for b\ and b 2 
(i.e., of the same order as the expected values for c\ and c 2 ). 

Since i t is not easy to disting uish the Zeeman signal in the V spectra of individual lines from the 



noise level, ICrutcher et al.l (J1999I ) produced an averaged sum of the N = 1 — >• hyperfine lines with 
strong Zeeman coefficients to display the Zeeman fit for their detection. In the case of the N = 2 — > 1 
transition, not all hyperfine lines have the same sign for their Zeeman coefficients and therefore it 
is only possible to average together the lines with same coefficient signs. The top panel of Figure 
[7] shows the weighted average Stokes I for lines 1, 2 and 4 that have negative Zeeman coefficients, 
with the weights being the relative sensitivity to the magnetic field, i.e., |Zj x R.Ij| given in Table [3j 
All the hyperfine lines are centered on the Vlsr of the source, i.e., -3 km s _1 and in order to remove 
the contamination of line 3 from line 2, the Gaussian fit for line 3 was subtracted from the blended 
spectrum. A similar procedure was performed for the averaged V spectrum of the aforementioned 
lines, which is displayed in the lower panel of Figure [7] with the instrumental effects removed. The 
bold line is the Zeeman fit to the average V data, expressed by V = Z{C\dIi/du) + Z{C 2 dI 2 jdv\ 
where Z is the weighted average of the negative Zeeman coefficients, with the weights being |R.Ij|. 
The fits for C\ and C 2 are C\ = —0.15 ± 2.25 mG and C 2 = —3.57 ± 3.22 mG, which are consistent 
with the values obtained for c\ and c 2 from the simultaneous Zeeman fit to individual V spectra 
mentioned above. 

The Zeeman Stokes V signal due to a field of a few hundred /iGauss is small (i.e., a line 
broadening of a few hundred Hz), and the noise level in the spectra needs to be sufficiently low 
in order to obtain at least a 2a detection. We developed simulations for the observed CN lines 
with different noise levels and magnetic field st r ength s, to estimate the remaining integration time 



required to obtain a detection. ICrutcher et al.l (|1999l ) obtained a value of B\ os ~ 0.75 mG for one 



of the velocity components of DR21(OH) by observing the N = 1 — > transition of CN, which is 
associated with a critical density n c ~ 10 5 cm~ 3 . Assuming that the magnetic field strength varies 
with n 1 ' 2 , where n is the the average gas density of the observed region, we should expect that the 
magnetic field we probe with our observations will be somewhat stronger. Since n c ~ 10 6 cm -3 at 
the N = 2 — > 1 transition, i?i os should be at most a few times stronger than the value obtained 
with the N = 1 — > transition. Our simulations show that for a 2 mG field, we need to have 
our noise level down to ~ 5 mK to be able to get at least a 2<r detection. The RMS noise in our 



10 



data is currently ~ 16 mK, implying that we need more observations to obtain a credible detection. 
Assuming an average system temperature of 400 K, with a bandwidth of 61 kHz, we will require 
about 48 more hour s of on-source integra tion time. This time estimate is comparable with the 



observing time that ICrutcher et al.l (|1999l ) spent on OMCln and DR210H to obtain a Zeeman 



detection in the CN (N = 1 — > 0) transition. 



5. Summary 

We recently designed and successfully commissioned a Four-Stokes-Parameter spectral line Po- 
larimeter (FSPPol) at the CSO in November 2008. The simple design of FSPPol does not contain 
any mirrors or grids to redirect or split the radiation beam, and the instrument is conveniently 
mounted in the elevation tube between the tertiary mirror behind the telescope dish and the Nas- 
myth focus where the heterodyne receiver is located. FSPPol transmits the beam to the receiver 
through half-wave and quarter-wave plates that are presently optimized for observations at 226 
GHz. 

We used FSPPol for linear and circular polarization measurements in the spectral lines of 
interstellar molecules during the months of July, September and October 2009. We measured a 
linear polarization level of ~ 1% to 2.5% due to the Goldreich-Kylafis effect in the spectral line wings 
of 12 CQ (J = 2 — > 1) in Orion KL/IRc2, and our results are consistent with previous observations 



(jGirart et al.ll2004l ). We also started Zeeman observations on the N = 2 — )■ 1 transition of CN in 
DR21(OH) for the first time. At this point we have obtained about 10 hours of on-source integration 
time, and our preliminary data analysis shows that although we have not detected a Zeeman signal 
in the CN (N = 2 — > 1) lines, the overall contribution from the instrumental effects in the Stokes 
V spectrum is low. Further observations of CN (N = 2 — > 1) are ongoing for this source and other 
star forming regions. 

The authors thank B. Dalrymple, S. H. Chen and T. Officer for their help in the develop- 
ment of FSPPoL, H. Shinnaga for helping with the Zeeman calculations and M. Azimlu and S. 
Chitsazzadeh for assistance during observations. M. H.'s research is funded through the NSERC 
Discovery Grant, Canada Research Chair, Canada Foundation for Innovation, Ontario Innovation 
Trust, and Western's Academic Development Fund programs. The Caltech Submillimeter Observa- 
tory is funded through the NSF grants AST 08-38261 and AST 05-40882 to the California Institute 
of Technology. 
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CSO elevation 
tube entrance 



Fig. 1. — The schematic diagram of FSPPol installed in the elevation tube at the CSO. The 
polarizing grid, used for testing, at the entrance of FSPPol is also shown. 
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Fig. 2. — The schematic diagram of the CSO telescope, also showing the position of FSPPol at 
1.04 m from the tertiary mirror towards Mirror 4. The telescope has an alt-azimuth mount and 
the receiver is located at one of the two Nasmyth foci. The virtual focus of the tertiary mirror is 
located 1.07 m behind the tertiary, and the beam waists (W) at the virtual focus and at FSPPol 
are 35.9 mm and 39.5 mm, respectively. 
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Fig. 3. — The graph of the variation of the power measured from a cold load with the CSO receiver 
against different orientations of the HWP with respect to the receiver polarization axis. The R x 
power is in arbitrary units. 
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e = (Y+e + 7t/2)/2 
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Fig. 4. — A schematic diagram for our linear polarization measurements. The equatorial coordinate 
system shows North and East on the sky. The receiver's polarization axis (R x ) projected on the sky 
in the equatorial system, lies along the East— West sky axis, 9 is the angle between R x and the slow 
axis of the HWP (denoted by "slow"), 7 is the parallactic angle of the source, and 9 the angle at 
which we want to measure the polarization state in the frame of the source. Accordingly, we have 
9 = ( 7 + 0' + 90°)/2. 



17 




Stokes Parameters - Sky Co 




Stokes Parameters - Sky Ct 




20 3C 



Fig. 5. — Linear polarization measurements of 12 CO (J = 2 — y 1) in Orion KL/IRc2; (top) Stokes 
/, corrected for the telescope beam efficiency; (middle) The polarization level p and (bottom) angles 
PA across the spectral line. 
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Fig. 6. — Stokes I spectrum of the nine strong hyperfine components of CN (N = 2 
DR21(OH). Line 5 is actually three blended lines that appear as one. 
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Fig. 7. — (top) The averaged Stokes I spectrum for the CN (N = 2—7-1) hyperfine lines 1, 2 & 4 
in DR21(OH) that have negative Zeeman coefficients. The solid lines are the Gaussian fits to the 
two velocity components of Stokes / spectrum, (bottom) The averaged Stokes V spectrum for the 
same hyperfine lines. The solid line is the fit for the Zeeman expression to the Stokes V spectrum, 
with instrumental effects removed. 
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Table 1 . Properties of the wave plates 



Property 



Material 

Indices of refraction: 

Ordinary 

Extraordinary 

Thickness 

Anti-reflection coating: 

Material 



single crystal Quartz 

2.106 ±0.006 a 

2.154 ± 0.007 a 

13.82 mm (HWP), 6.90 mm (QWP) 

High density Polypropylene (HDPP) 



a Birch et al. 1994, IEEE: Transactions on Microwave Theory 
and Techniques, 42, 6, 956 
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Table 2. Polarization measurements on CO J = 2 — > 1 in Orion KL/IRc2 



_^ Pj%) PA (deg) 

Offsets (arcsec) a Blue b wing Line center Red c wing Blue b wing Red c wing 

(20", 20") 2.1 ± 0.3 0.35 ± 0.04 1.8 ± 0.2 77.1 ± 4.2 81.2 ± 3.4 

(20", -20") 1.6 ± 0.2 0.29 ± 0.04 2.4 ± 0.2 97.8 ± 3.3 88.1 ± 2.9 

(-20", -20") 1.1 ± 0.3 0.29 ± 0.07 1.5 ± 0.4 154.8 ± 7.1 155.6 ± 8.1 

a Offsets are with respect to IRc2 (a = 05 h 35 m 14.5 s , 8 = -05°22'30.4"; J2000.0). 
b p and PA are averaged over the velocity range of -4 to 6 km s _1 (see Figure(5]). 
c p and PA are averaged over the velocity range of 13 to 23 km s _1 . 
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Table 3. Hyperfine components of CN (N = 2 — y 1) 



Line 


(N,J,F)^(N',J',F') 


v (MHz) 


Z (Hz/^G) 


R.I b 


|ZxR.I| c 


1 


(2,3/2,3/2)^(1,1/2,3/2) 


226632.19 


-0.72241 


0.59259 


0.42809 


2 


(2,3/2,5/2)^(1,1/2,3/2) 


226659.58 


-0.70995 


2.0 


1.41991 


3 


(2,3/2,1/2)^(1,-1/2,1/2) 


226663.70 


0.62277 


0.59259 


0.36905 


4 


(2,3/2,3/2)^(1,1/2,1/2) 


226679.38 


-1.18326 


0.74074 


0.87649 


5 


(2,5/2,5/2)^(1,3/2,3/2) 


226874.17 


0.70995 


2.016 


1.43127 


5 


(2,5/2,7/2)^(1,3/2,5/2) 


226874.75 


0.40035 


3.200 


1.28112 


5 


(2,5/2,3/2)^(1,3/2,1/2) 


226875.90 


1.18326 


1.200 


1.41991 


6 


(2,5/2,3/2)^(1,3/2,3/2) 


226887.35 


1.46973 


0.3840 


0.56438 


7 


(2,5/2,5/2)^(1,3/2,5/2) 


226892.12 


1.05692 


0.3840 


0.40586 



a Shinnaga et al. (2010), in preperation. 
b Relative intensity. 
c Relative sensitivity to B\ os . 



